DecodeEPSG Function

public function DecodeEPSG(epsg) result(gridmapping)

get EPSG code as input and set coordinate reference system

Arguments

Type IntentOptional Attributes Name
integer(kind=short), intent(in) :: epsg

Return Value type(CRS)


Variables

Type Visibility Attributes Name Initial
real(kind=float), public :: centralMeridian

Source Code

FUNCTION DecodeEPSG &
!
(epsg) &
!
RESULT (gridmapping)


IMPLICIT NONE

!Arguments with intent in:
INTEGER (KIND = short), INTENT(IN) :: epsg


!Local declarations:
TYPE (CRS) :: gridmapping
REAL (KIND = float)   :: centralMeridian

!----------------------end of declarations-------------------------------------
gridmapping % epsg = epsg
SELECT CASE (epsg)

CASE(4326)
  CALL SetCRS (GEODETIC, WGS84, gridmapping)
  CALL SetGeodeticParameters (gridmapping, prime_meridian = 0.0)

CASE(4230)
  CALL SetCRS (GEODETIC, ED50, gridmapping)
  CALL SetGeodeticParameters (gridmapping, prime_meridian = 0.0)

CASE(3003) !Gauss Boaga West
  CALL SetCRS (TM, ROME40, gridmapping)
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = 9. * degToRad, &
              falseE = 1500000., falseN = 0., k = 0.9996)

CASE(3004) !Gauss Boaga East
  CALL SetCRS (TM, ROME40, gridmapping)
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = 15. * degToRad, &
              falseE = 2520000., falseN = 0., k = 0.9996)

CASE(21781) !CH1903 Swiss topo
  CALL SetCRS (SOC, CH1903, gridmapping)
  CALL SetSwissParameters &
          (gridmapping, latc = 0.819474, lonc = 0.129845, &
           azimuth = 1.570796, falseE = 600000., falseN = 200000., k = 1.)

CASE(23028) !ED50 / UTM zone 28N
  CALL SetCRS (TM, ED50, gridmapping)
  centralMeridian = (6 * 28 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(23029) !ED50 / UTM zone 29N
  CALL SetCRS (TM, ED50, gridmapping)
  centralMeridian = (6 * 29 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(23030) !ED50 / UTM zone 30N
  CALL SetCRS (TM, ED50, gridmapping)
  centralMeridian = (6 * 30 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(23031) !ED50 / UTM zone 31N
  CALL SetCRS (TM, ED50, gridmapping)
  centralMeridian = (6 * 31 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(23032) !ED50 / UTM zone 32N
  CALL SetCRS (TM, ED50, gridmapping)
  centralMeridian = (6 * 32 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(23033) !ED50 / UTM zone 33N
  CALL SetCRS (TM, ED50, gridmapping)
  centralMeridian = (6 * 33 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(23034) !ED50 / UTM zone 34N
  CALL SetCRS (TM, ED50, gridmapping)
  centralMeridian = (6 * 34 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(23035) !ED50 / UTM zone 35N
  CALL SetCRS (TM, ED50, gridmapping)
  centralMeridian = (6 * 35 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(23036) !ED50 / UTM zone 36N
  CALL SetCRS (TM, ED50, gridmapping)
  centralMeridian = (6 * 36 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(23037) !ED50 / UTM zone 37N
  CALL SetCRS (TM, ED50, gridmapping)
  centralMeridian = (6 * 37 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(23038) !ED50 / UTM zone 38N
  CALL SetCRS (TM, ED50, gridmapping)
  centralMeridian = (6 * 38 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32601) !WGS 84 / UTM zone 1N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 1 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32602) !WGS 84 / UTM zone 2N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 2 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32603) !WGS 84 / UTM zone 3N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 3 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32604) !WGS 84 / UTM zone 4N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 4 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32605) !WGS 84 / UTM zone 5N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 5 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32606) !WGS 84 / UTM zone 6N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 6 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32607) !WGS 84 / UTM zone 7N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 7 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32608) !WGS 84 / UTM zone 8N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 8 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32609) !WGS 84 / UTM zone 9N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 9 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32610) !WGS 84 / UTM zone 10N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 10 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32611) !WGS 84 / UTM zone 11N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 11 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32612) !WGS 84 / UTM zone 12N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 12 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32613) !WGS 84 / UTM zone 13N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 13 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32614) !WGS 84 / UTM zone 14N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 14 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32615) !WGS 84 / UTM zone 15N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 15 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32616) !WGS 84 / UTM zone 16N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 16 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32617) !WGS 84 / UTM zone 17N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 17 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32618) !WGS 84 / UTM zone 18N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 18 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32619) !WGS 84 / UTM zone 19N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 19 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32620) !WGS 84 / UTM zone 20N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 20 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32621) !WGS 84 / UTM zone 21N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 21 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32622) !WGS 84 / UTM zone 22N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 22 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32623) !WGS 84 / UTM zone 23N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 23 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32624) !WGS 84 / UTM zone 24N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 24 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32625) !WGS 84 / UTM zone 25N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 25 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32626) !WGS 84 / UTM zone 26N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 26 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32627) !WGS 84 / UTM zone 27N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 27 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32628) !WGS 84 / UTM zone 28N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 28 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32629) !WGS 84 / UTM zone 29N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 29 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32630) !WGS 84 / UTM zone 30N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 30 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32631) !WGS 84 / UTM zone 31N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 31 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32632) !WGS 84 / UTM zone 32N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 32 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32633) !WGS 84 / UTM zone 33N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 33 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32634) !WGS 84 / UTM zone 34N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 34 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32635) !WGS 84 / UTM zone 35N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 35 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32636) !WGS 84 / UTM zone 36N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 36 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32637) !WGS 84 / UTM zone 37N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 37 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32638) !WGS 84 / UTM zone 38N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 38 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32639) !WGS 84 / UTM zone 39N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 39 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32640) !WGS 84 / UTM zone 40N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 40 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32641) !WGS 84 / UTM zone 41N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 41 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32642) !WGS 84 / UTM zone 42N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 42 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32643) !WGS 84 / UTM zone 43N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 43 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32644) !WGS 84 / UTM zone 44N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 44 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32645) !WGS 84 / UTM zone 45N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 45 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32646) !WGS 84 / UTM zone 46N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 46 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32647) !WGS 84 / UTM zone 47N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 47 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32648) !WGS 84 / UTM zone 48N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 48 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32649) !WGS 84 / UTM zone 49N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 49 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32650) !WGS 84 / UTM zone 50N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 50 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32651) !WGS 84 / UTM zone 51N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 51 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32652) !WGS 84 / UTM zone 52N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 52 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32653) !WGS 84 / UTM zone 53N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 53 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32654) !WGS 84 / UTM zone 54N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 54 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32655) !WGS 84 / UTM zone 55N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 55 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32656) !WGS 84 / UTM zone 56N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 56 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32657) !WGS 84 / UTM zone 57N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 57 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32658) !WGS 84 / UTM zone 58N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 58 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32659) !WGS 84 / UTM zone 59N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 59 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32660) !WGS 84 / UTM zone 60N
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 60 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 0., k = 0.9996)

CASE(32701) !WGS 84 / UTM zone 1S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 1 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32702) !WGS 84 / UTM zone 2S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 2 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32703) !WGS 84 / UTM zone 3S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 3 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32704) !WGS 84 / UTM zone 4S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 4 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32705) !WGS 84 / UTM zone 5S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 5 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32706) !WGS 84 / UTM zone 6S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 6 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32707) !WGS 84 / UTM zone 7S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 7 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32708) !WGS 84 / UTM zone 8S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 8 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32709) !WGS 84 / UTM zone 9S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 9 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32710) !WGS 84 / UTM zone 10S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 10 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32711) !WGS 84 / UTM zone 11S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 11 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32712) !WGS 84 / UTM zone 12S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 12 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32713) !WGS 84 / UTM zone 13S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 13 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32714) !WGS 84 / UTM zone 14S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 14 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32715) !WGS 84 / UTM zone 15S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 15 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32716) !WGS 84 / UTM zone 16S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 16 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32717) !WGS 84 / UTM zone 17S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 17 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32718) !WGS 84 / UTM zone 18S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 18 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32719) !WGS 84 / UTM zone 19S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 19 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32720) !WGS 84 / UTM zone 20S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 20 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32721) !WGS 84 / UTM zone 21S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 21 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32722) !WGS 84 / UTM zone 22S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 22 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32723) !WGS 84 / UTM zone 23S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 23 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32724) !WGS 84 / UTM zone 24S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 24 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32725) !WGS 84 / UTM zone 25S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 25 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32726) !WGS 84 / UTM zone 26S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 26 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32727) !WGS 84 / UTM zone 27S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 27 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32728) !WGS 84 / UTM zone 28S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 28 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32729) !WGS 84 / UTM zone 29S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 29 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32730) !WGS 84 / UTM zone 30S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 30 + 177) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32731) !WGS 84 / UTM zone 31S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 31 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32732) !WGS 84 / UTM zone 32S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 32 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32733) !WGS 84 / UTM zone 33S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 33 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32734) !WGS 84 / UTM zone 34S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 34 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32735) !WGS 84 / UTM zone 35S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 35 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32736) !WGS 84 / UTM zone 36S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 36 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32737) !WGS 84 / UTM zone 37S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 37 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32738) !WGS 84 / UTM zone 38S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 38 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32739) !WGS 84 / UTM zone 39S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 39 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32740) !WGS 84 / UTM zone 40S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 40 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32741) !WGS 84 / UTM zone 41S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 41 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32742) !WGS 84 / UTM zone 42S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 42 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32743) !WGS 84 / UTM zone 43S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 43 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32744) !WGS 84 / UTM zone 44S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 44 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32745) !WGS 84 / UTM zone 45S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 45 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32746) !WGS 84 / UTM zone 46S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 46 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32747) !WGS 84 / UTM zone 47S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 47 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32749) !WGS 84 / UTM zone 49S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 49 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32750) !WGS 84 / UTM zone 50S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 50 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32751) !WGS 84 / UTM zone 51S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 51 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32752) !WGS 84 / UTM zone 52S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 52 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32753) !WGS 84 / UTM zone 53S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 53 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32754) !WGS 84 / UTM zone 54S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 54 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32755) !WGS 84 / UTM zone 55S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 55 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32756) !WGS 84 / UTM zone 56S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 56 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32757) !WGS 84 / UTM zone 57S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 57 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32758) !WGS 84 / UTM zone 58S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 58 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32759) !WGS 84 / UTM zone 59S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 59 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

CASE(32760) !WGS 84 / UTM zone 60S
  CALL SetCRS (TM, WGS84, gridmapping)
  centralMeridian = (6 * 60 - 183) * degToRad
  CALL SetTransverseMercatorParameters &
             (gridmapping, lat0 = 0., centM = centralMeridian, &
              falseE = 500000., falseN = 10000000., k = 0.9996)

END SELECT

RETURN
END FUNCTION DecodeEPSG